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Abstract 



Results are presented from a search for long-lived neutralinos decaying into a 
photon and an invisible particle, a signature associated with gauge-mediated super- 
symmetry breaking in supersymmetric models. The analysis is based on a 4.9 fb -1 
sample of proton-proton collisions at a/s = 7 TeV, collected with the CMS detector 
at the LHC The missing transverse energy and the time of arrival of the photon at 
the electromagnetic calorimeter are used to search for an excess of events over the 
expected background. No significant excess is observed, and lower limits at the 95% 
confidence level are obtained on the mass of the lightest neutralino, m (x^) > 220 GeV 
(for ct < 500 mm), as well as on the proper decay length of the lightest neutralino, 
ct > 6000 mm (for m{jft) < 150 GeV). 
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1 Introduction 

New, heavy particles with long lifetimes are predicted in many models of physics beyond 
the standard model (SM), such as hidden valley scenarios U or supersymmetry (SUSY) with 
gauge-mediated supersymmetry breaking (GMSB) J2J. Under the assumption of R-parity con- 
servation [3j, strongly-interacting supersymmetric particles are pair-produced at the Large 
Hadron Collider (LHC). The decay chain may include one or more quarks and gluons, as well 
as the lightest supersymmetric particle (LSP), which escapes detection, giving rise to a momen- 
tum imbalance in the transverse plane. A GMSB benchmark scenario, commonly described as 
'Snowmass Points and Slopes 8' (SPS8) (U is used as the reference in this search. In this sce- 
nario, the neutralino is the next-to-lightest supersymmetric particle. It decays to a photon (or a 
Z boson) and a gravitino, which is the LSP 0- 

In each signal event, two neutralinos are produced, and therefore one or two photons in the 
majority of events. The search criteria require only one identified photon in order to be sen- 
sitive to scenarios with a large branching fraction for the neutralino decay to a Z boson and a 
gravitino. This strategy also provides sensitivity to models that have only one long-lived par- 
ticle that can decay to a photon. For a long-lived neutralino, the photon from the — > 7 G 
decay is produced at the decay vertex, at some distance from the beam line, and reaches the 
detector at a later time than the prompt, relativistic particles produced at the interaction point. 
In addition, the geometric shape of the energy deposit produced by such photons is typically 
different from that of a prompt photon. The time of arrival of the photon at the detector and 
the missing transverse energy are used to discriminate signal from background. 

A search for a long-lived neutralino, decaying to a photon and a gravitino, is performed with 
a novel technique using the excellent time measurement with the electromagnetic calorimeter. 
Previous searches for long-lived neutralinos have been performed by the CMS collaboration (6|, 
using the impact parameter of converted photons relative to the beam collision point, and by 
the CDF [7] collaboration. Other searches with prompt photons by the ATLAS |8] and DO [9] 
collaborations place lower limits on the mass of the at 280 GeV and 175 GeV, respectively, in 
the SPS8 scenario. 

2 Detector and data samples 

A detailed description of the Compact Muon Solenoid (CMS) detector can be found else- 
where flT0| . The detector's central feature is a superconducting solenoid providing a 3.8 T 
axial magnetic field along the beam direction. Charged particle trajectories are measured by 
a silicon pixel and strip tracker system with full azimuthal coverage within \t]\ < 2.5; the 
pseudo-rapidity rj is defined as rj = — In [tan (9/2)], with 8 being the polar angle with respect 
to the counterclockwise beam direction. A lead-tungstate crystal electromagnetic calorime- 
ter (ECAL) and a brass/ scintillator hadron calorimeter (HCAL) surround the tracker volume. 
The ECAL is a high-granularity device constructed using PbWC>4 crystals. The barrel region 
consists of 61 200 crystals with a frontal area of approximately 2.2 x 2.2 cm 2 corresponding to 
roughly 0.0174 x 0.0174 in rj-cp space. Each of the two endcap sections consist of 3662 crystals 
with a frontal area of 2.68 x 2.68 cm 2 . A typical shower spans approximately 10 crystals with 
energy deposits above the threshold. Muons are measured in gas-ionization detectors embed- 
ded in the steel return yoke of the magnet. The detector is nearly hermetic, allowing reliable 
measurement of transverse momentum imbalance to be performed. The time of arrival of elec- 
tromagnetic particles can be measured to excellent precision using the CMS ECAL fllT|. The 
time reconstruction method is described in more detail in Section I37T1 
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3 Analysis technique 



The analysis is performed on the proton-proton collision data at a center-of-mass-energy of 
7 TeV recorded by the CMS detector at the LHC, corresponding to an integrated luminosity of 
4.9 ± 0.1 fb _1 . Events with at least one high transverse momentum (px) isolated photon in the 
barrel region and at least three jets in the final state are selected in this analysis. The data were 
recorded using the CMS two-level trigger system. Several trigger selections have been used 
due to the increasing instantaneous luminosity during the data taking. The first 0.20 fb -1 of 
data were collected with a trigger requiring at least one isolated photon with pj > 75 GeV. For 
the second 3.8 fb -1 , the pj threshold was increased to 90 GeV. In the remaining 0.89 fb _1 , the 
trigger selection required at least one isolated photon with pj > 90 GeV in the barrel region 
and at least three jets with pj greater than 25 GeV. All offline selection requirements are chosen 
to be more restrictive than the trigger selection. 

Signal and background events are generated using Monte Carlo (MC) packages pythia 6.4.22 
ff!2| or MadGraph 5 E3| with the CTEQ6L1 HU parton distribution functions (PDFs). The re- 
sponse of the CMS detector is simulated using the Geant4 package 1 15 1 . Decays of secondary 
t leptons, coming from W and Z productions, are handled with TAUOLA I1T61 . The SUSY GMSB 
signal production follows the SPS8 proposal, where the free parameters are the SUSY breaking 
scale (A) and the average proper decay length (ct) of the neutralino. The mass explored is in 
the range of 140 to 260 GeV (corresponding to A values from 100 to 180 TeV), with proper decay 
lengths ranging from ct = 1 mm to 6000 mm. These free parameters are varied to cover the 
range of experimental phase space allowed by inner radius of the barrel section of the ECAL 
(1.29 m). 

There is a non-negligible probability that several collisions may occur in a single bunch crossing 
due to the high instantaneous luminosities at the LHC. The presence of multiple interaction 
vertices in an event (pile-up) affects the resolution of the transverse momentum measurement 
and the performance of photon isolation requirements. To account for the effects of pile-up, 
simulated events are re-weighted so that the distribution of the number of interaction vertices 
matches that in the data. 

3 Analysis technique 

This section, outlining the analysis technique, starts with a description of the physics object 
reconstruction followed by a brief explanation of the event selection criteria. Finally, the def- 
initions of the key discriminating variables related to the ECAL cluster shape and the time of 
impact of the photon on the surface of the ECAL are discussed. The signal and background 
yields are determined with a binned maximum likelihood fit to the two-dimensional distribu- 
tion in these variables. 

3.1 Object reconstruction 

Photons are reconstructed by identifying energy deposits in the ECAL using the method ex- 
plained in Ref. ffl7|. Photons that are found to have converted into an electron-positron pair 
in the detector material are not used in the analysis. Electron candidates are reconstructed 
starting from a cluster of energy deposits in the ECAL which is then matched to the momen- 
tum associated with a track in the silicon tracker. Electron candidates are required to have 
\t] | < 1.44 or 1.56 < \t]\ < 2.5 to avoid the region of transition between the barrel and 
endcap sections. Photons are required to be spatially separated from electrons by at least 
AR = (A?/) 2 + (A(p) 2 = 0.25, where Arj and A(p are, respectively, differences between the 
photon and the electron directions in pseudorapidity and azimuthal angle. 



3.1 Object reconstruction 
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Jets are reconstructed from objects identified using the Particle-Flow (PF) algorithm [18 j with 
anti-fcx clustering [19] and a distance parameter of 0.5. In this analysis, the missing transverse 
energy (£x) is defined as the magnitude of the vector sum of the transverse momentum of all 
particles identified in the PF algorithm in the event excluding muons. 

The time of impact, T raw , for the photon on the surface of the ECAL is the weighted time of 
impact measured in the crystals within the cluster associated with a photon candidate. An 
event-by-event correction (T promp t) is applied to T raw to account for possible biases due to the 
jitter in the trigger system, and to the imperfect knowledge of the time of the interaction within 
the bunch crossing. This correction is computed using the time of impact of all crystals in 
the event, excluding those belonging to the two most energetic photon candidates, which are 
typically due to prompt jets, low-energy prompt photons, and photons from 7T° and rj decays. 
The new calibrated ECAL timing is defined as T ca ]ib — T raw — Tp rom pt- With this definition, 
a particle produced at the interaction point has a time of arrival of zero, whereas a delayed 
photon has a non-zero T ca iib- The distributions in data for T raw and T ca iib, after the nominal 
selection, are shown in Fig. [T] The width of the main, Gaussian, component of T ca iit, is slightly 
smaller than that of T raw , while there is some increase in the tails. For the dominant background 
processes, the tails are taken into account by using control samples in data, as described in 
Section|4] In the determination of the yield, the distribution of T ca m, in simulated signal events is 
used as a template for the signal contribution. This distribution is narrower in simulation than 
in the data, because the uncertainties in the time inter-calibration constants are not emulated. 
A convolution with a Gaussian, whose parameters vary as a function of the photon energy, is 
performed to reproduce the T ca m, resolution observed in data. 



ECAL Timing [ns] 

Figure 1: The ECAL timing distribution for data, before and after calibration, overlaid with the 
results of the Gaussian fits. 

One of the distinctive features of a photon is the shape of the energy deposits it leaves in the 
ECAL. Prompt photons have a roughly circular projected energy deposit on the ECAL sur- 
face, while the energy deposits from jets typically have a larger width along the rj direction. 
Non-prompt photons are expected to have an elliptical shape along an arbitrary direction, as 
illustrated in Fig. |2j therefore the width of the energy deposit along the fj direction is not opti- 
mal for the discrimination of jets. In this search, the shape of the energy deposit is characterized 
by the minor axis (SMinor) of its projection on the internal ECAL surface. The axis SMm r is com- 
puted using the geometrical properties of the distribution of the energy deposit, and is defined 



CMS Data 



\fs = 7 TeV 



c 30000 




4 



3 Analysis technique 



as 



(1) 



where S^, S^, and are the second moments of the spatial distribution of the energy deposit 
in the ECAL in tj-(p coordinates. A large fraction of QCD multijet events is rejected by applying 
requirements on SMinor as illustrated in Fig. |3} where the normalized distributions of SMinor for 
simulated signal and QCD multijet background events are shown. 

3.2 Event selection 

Events must have a primary vertex with at least four associated tracks and a position less 
than 2 cm from the center of the CMS detector in the direction transverse to the beam and 
24 cm in the direction along the beam. Events are also required to have at least three jets with 
pj > 35 GeV and spatially separated from photons by at least AR = 0.5. 

Photon candidates are required to have pj > 100 GeV and \tj\ < 1.4 and to be isolated in the 
HCAL, the ECAL, and the tracker. An absolute isolation parameter is defined as the scalar sum 
of the transverse energies of tracks or calorimeter deposits in a cone of aperture 0.3 around 
the photon direction, excluding the contribution from the photon itself. A relative isolation 
parameter is defined as the ratio of the absolute isolation and the photon Ej. In the tracker, the 
relative isolation is required to be less than 0.1. In the ECAL and the HCAL, the relative isola- 
tion is required to be less than 0.05 and the absolute isolation less than 2.4 GeV. Thresholds on 
both absolute and relative isolation are set in the ECAL and HCAL to avoid imposing require- 
ments that are more restrictive than the noise level. The energy deposit by a photon candidate 
is required to have 0.15 < SMinor < 0-30. This requirement is optimized to select candidates 
that are more likely to be real photons. 

The signal efficiencies for selecting one photon and at least three jets are summarized in Table [l] 
for proper decay lengths between 1 mm and 6000 mm and for A between 100 TeV and 180 TeV. 
The efficiency drops by a factor of two between ct = 1 mm and 6000 mm, since, with increasing 
decay time, the probability of the to decay outside the detector is enhanced. 
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Figure 2: The distribution of energy deposition in the ECAL crystals for a prompt (left) and a 
non-prompt (right) photon. Each rectangle represents an ECAL crystal and has a size that is 
proportional to the energy deposited in that crystal. 
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Figure 3: Normalized distribution of SMinor for simulated signal and QCD multijet events. 



Table 1: Selection efficiency in percent. The reported uncertainties include the contributions of 
systematic effects, for various signal samples. 



A (TeV) 


M^o (GeV) 


ct = 1 mm 


ct = 250 mm 


ct = 2000 mm 


ct = 6000 mm 


100 


140 


18.7 ± 0.3 


18.4 ± 0.2 


8.4 ± 0.2 


3.3 ± 0.1 


120 


170 


24.9 ± 0.4 


24.6 ± 0.3 


15.1 ± 0.4 


6.6 ± 0.1 


140 


200 


30.4 ± 0.3 


31.3 ± 0.3 


22.2 ± 0.4 


11.4 ± 0.3 


160 


230 


35.5 ± 0.3 


36.1 ± 0.6 


29.4 ± 0.4 


17.0 ± 0.4 


180 


260 


40.1 ± 0.7 


38.0 ± 0.5 


36.0 ± 0.5 


22.2 ± 0.4 



4 Background estimation 



The primary sources of background in the analysis are QCD multijet events and 7 + jets events, 
which together make up 99% of the sample. Improper reconstruction of jets can give rise to fake 
missing transverse energy, while photons produced in the decays of hadrons (mostly energy 
7T° and rj) can sometimes pass the isolation criteria. 

A large fraction of 7 + jets events, characterized by a smaller jet multiplicity compared to signal, 
are rejected by requiring at least three jets in the event. The residual contribution of these 
backgrounds is estimated from the data. 

In addition, there are other (non-QCD) processes with genuine fij, largely comprised of W/ Z + 
7 + jets and tt events, where the W boson decays into a lepton and a neutrino. There is also a 
small contribution from Drell-Yan processes. These events make up less than 1% of the total 
sample but are taken into account since they can play a role in the tails of the £t distribution 
where signal is expected. Simulated events are used to estimate the contribution of these pro- 
cesses. 

Finally, additional backgrounds from events not originating from proton-proton collisions, in- 
cluding cosmic rays and beam-halo muons, are also expected. The contribution of these events 
is reduced to negligible levels by requiring T ca \fo of the most energetic photon candidate to be 
greater than —2 ns, and the event to have an identified primary vertex and at least three jets. 

Because of the difficulty of accurately predicting cross sections and jet multiplicities for multijet 
and 7 + jets processes, their contribution is estimated with methods based on the data. The 
QCD multijet control sample is obtained by selecting events with at least three jets and a photon 
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5 Systematic uncertainties 



candidate passing a less stringent identification requirement but failing the nominal photon 
selection criteria. The 7 + jets control sample consists of events with one photon which satisfies 
the nominal selection and with the highest-px jet (leading jet) back-to-back with respect to the 
photon, in the transverse plane. The ratio of the transverse momenta of the leading jet to 
that of the photon is required to be between 0.6 and 1.4, while for the subleading jet the ratio 
is required to be less than 0.2. The contribution of non-QCD and signal events to these two 
control samples is estimated to be, respectively, 1% and less than 0.01%. 

To estimate the number of background and signal events in data, a maximum likelihood fit is 
performed to the distributions of and T ca \^, which are uncorrected. Binned shape templates 
are derived from simulated events for signal and non-QCD backgrounds. Templates for QCD 
multijet and 7 + jets are derived from data control samples as described earlier. The relative 
normalization of the QCD multijet and 7 + jets components is fixed to 67% and 33%, respec- 
tively, based on studies with simulated events. The normalization of the non-QCD templates 
are fixed in the fit according to the measured cross sections (statistical uncertainties in the cross 
sections are less than 3%) and the integrated luminosity of the data sample. Studies have been 
performed with pseudo-experiments to confirm the stability of the fit and to verify that the fit 
results are unbiased. The measured signal and background yields in data, obtained with the 
likelihood fit, are summarized in Table [2] The one-dimensional projections of £x and T ca \i\, for 
the data and expected backgrounds, as determined from the fit, are illustrated in Fig. |4j No 
excess of events is observed beyond the SM backgrounds and the fitted signal yield is compat- 
ible with zero. It should be noted that the discriminating power of individual variables is not 
apparent in these projections because the largest sensitivity to signal is in the region with both 
large £x and large T ca \sb- The improved background discrimination is visible in Fig.[5]where the 
one-dimensional projection of fij for events with T ca \Q, > 0.5 ns is illustrated. 

Table 2: The measured signal and background yields determined with the maximum likeli- 
hood fit to the data. The relative composition of QCD multijet and 7 + jets backgrounds have 
been normalized to 67% and 33% with respect to each other. The expected signal yields are, 
respectively, 211 events for GMSB(100,250) and 96 for GMSB( 100,2000) benchmark points. The 
reported uncertainties are statistical only and are determined in the fit. 



Events 


GMSB (100, 250) 


6±8 


GMSB (100, 2000) 


4±4 


QCD multijet and 7 + jets 


80900 ± 300 


tt + jets (fixed) 


73 


W — > ev + jets (fixed) 


116 


Drell-Yan + jets (fixed) 


67 


W/Z+ jets + 7 (fixed) 


215 


Total background 


81400 ± 300 


Data 


81382 



5 Systematic uncertainties 

Several sources of systematic uncertainty have been considered and their contributions are 
summarized in Table [3] The largest single contribution to the systematic uncertainties derives 
from the uncertainty in the modeling of the background shape. A bin-by-bin variation of the 
background shape template according to the Poisson uncertainty is used to determine the con- 
tribution of each type of background. An additional uncertainty is assessed for the QCD mul- 
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Figure 4: The one-dimensional projection for Ej (left) and for ECAL timing (right), after all 
selection requirements. The multijet and 7 + jets backgrounds are normalized to the yields 
from the fit. The rest of the backgrounds are fixed according to the integrated luminosity of the 
data. The GMSB(100,2000) benchmark point corresponds to A = 100 TeV, cr = 2000 mm and 
the GMSB(100,250) benchmark point corresponds to A = 100 TeV, ct = 250 mm. 



tijet and 7 + jets processes using simulated events, by comparing the shapes of Ej and 
for the control sample and and for a sample obtained with the nominal selection criteria. The 
difference observed in simulation is used to re-weight the shapes obtained in data control sam- 
ples. The dominant contribution is due to the difference in the Ej distributions. The small 
tails in the distribution of T ca ub are accounted for by using data control samples to derive the 
templates, rather than relying on simulation. The uncertainty in the relative fraction of QCD 
multijet and 7 + jets events is taken to be equal to observed difference of 10% between the 
number of events in the 7 + jets control sample in data and the expected number of events 
according to PYTHIA." 

The main contributions to the uncertainty in the signal shape modeling derive from the uncer- 
tainty in the #T resolution and the determination of T ca Hb. The contribution of the j£j resolution 
uncertainty is estimated by smearing the Ej distribution of simulated signal events. A sys- 
tematic uncertainty of 0.1 ns is assigned to the measurement of the time of impact T ca \fo. This 
uncertainty is determined using a sample of 7 + jets events by measuring the difference be- 
tween the average T ca iib values in data and simulation, as a function of the photon Ej. 

The uncertainty in the luminosity determination is 2.2% 1120 1 . The remaining sources of system- 
atic uncertainty affecting the signal acceptance are the following. The calorimeter response to 
different types of particles is not perfectly linear and hence corrections are made to properly 
map the measured jet energy deposition. The uncertainty on this correction is referred to as the 
uncertainty on the jet energy scale and varies as a function of position and transverse momen- 
tum of the jet. Similarly, the uncertainty on the photon energy scale in the barrel is estimated 
to be 1.0%, based on the final-state radiation measurement with Z bosons fl2"T| . Following the 
recommendations of the PDF4LHC group [22J, PDF and the strong coupling constant (a s ) vari- 
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6 Results 
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Figure 5: The one-dimensional projection for j^t after all selection requirements, for events 
with T ca \Q, > 0.5 ns. The multijet and 7 + jets backgrounds are normalized to the yields from 
the fit. The rest of the backgrounds are fixed according to the integrated luminosity of the data. 
The GMSB(100,2000) benchmark point corresponds to A = 100 TeV, ct = 2000 mm and the 
GMSB(100,250) benchmark point corresponds to A = 100 TeV, ct = 250 mm. 



ations of the MSTW2008 |23|, CTEQ6.6 EH and NNPDF2.0 [25 1 PDF sets are taken into account 
and their impact on the signal acceptance is estimated. 



6 Results 

The observed event yield in data is consistent with the SM background prediction, and upper 
limits are obtained on the production cross section of a long-lived neutralino in the context of 
the GMSB model, assuming BQg — > 7G) = 100%. Exclusion limits are computed with a mod- 
ified frequentist CL S method [26-28], using the asymptotic approximation for the test statistic 
as described in Ref. |29| . The background normalization and the corresponding uncertainty are 
taken from the fit to the data. The uncertainties in the shapes are taken into account by vertical 
interpolation of the templates. The shapes are interpolated quadratically for shifts below one 
standard deviation and linearly beyond. Log-normal multiplicative corrections are used for the 
normalization, the signal acceptance, and the integrated luminosity. Fig. [6] shows the observed 
and expected 95% confidence level (CL) upper limits on the cross section for GMSB production 
in terms of mass (left), and proper decay length (right). The signal cross section is computed 
at leading order precision and the theoretical uncertainty is evaluated by using the PDF4LHC 
recommendation for the PDF uncertainty. The one-dimensional limits are combined to provide 
exclusion limits in the mass and proper decay length plane of the x\ i R Fig-H 
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Table 3: Summary of the systematic uncertainties in the background and signal shapes, as well 
as in the signal acceptance x efficiency. The signal uncertainties are evaluated individually 
for every signal point, although only the maximum and minimum values associated with each 
source are quoted. 



Source 


Uncertainty (%) 


Background 


Shape 


10 


Normalization 


0.3 


Multijet/ 7 + jets fraction 


0.5 


Signal shape 


#r resolution 


0.2-2 


ECAL timing uncertainty 


1-5 


Signal acceptance x 


efficiency 


Photon energy scale 


0.5-3 


Jet energy scale 


0.02-0.05 


Jet energy resolution 


0.01-2 


PDF uncertainties 


0.1-2 
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170 GeV (right) in the SPS8 



7 Summary 

The CMS experiment has performed a search for long-lived particles produced in association 
with jets using LHC proton-proton collision data at a center-of-mass energy of 7 TeV corre- 
sponding to an integrated luminosity of 4.9 ± 0.1 fb -1 . A GMSB scenario with a long-lived 
neutralino decaying to a photon and a gravitino is used as the reference. The missing trans- 
verse energy and the timing information from the ECAL are used to search for an excess of 
events over the expected SM background prediction. A fit to the two-dimensional distribution 
in these variables yields no significant excess of events beyond the SM contributions, and up- 
per limits at 95% CL are obtained on the GMSB production cross section in the SPS8 model 
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7 Summary 
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SPS8 model of GMSB supersymmetry. 



of GMSB supersymmetry. In this scheme, we obtain an exclusion region as a function of both 
the neutralino mass and its proper decay length. The mass of the lightest neutralino is then 
restricted to values tn{j^) > 220 GeV (for neutralino proper decay length ct < 500 mm) at 95% 
CL, and the neutralino decay length ct must be greater than 6000 mm (for tn{]^) < 150 GeV). 
These limits are the most stringent for long-lived neutralinos. 
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